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SPIN AND RECOVERY CHARACTERISTICS OF 
THE NORTH AMERICAN X-15 AIRPLANE 
By William D. Grantham and Stanley H. Scher 

SUMMARY 


An analytical study has been made to aid in predicting the spin and 
recovery characteristics of the North American X-15 atmospheric-reentry 
research airplane. The approach used was to simulate first on a high- 
speed digital computer the spin entry, the ensuing spin, and the recovery 
therefrom obtained from dynamic-model tests at Low Reynolds number by 
using appropriate aerodynamic force and moment data. Then, high Reynolds 
number aerodynamic data were inserted into the computer and the resulting 
effects oh spin entry and spin characteristics were calculated. 

The results indicated that the airplane with the lower rudder on or 
off is not likely to enter a spin from trimmed gliding or level flight. 

A developed spin might be obtained if the airplane receives some violent 
disturbance which would tend to put it in a flight condition having rota- 
tion and a high angle of attack, approximating the attitude and motion of 
a model when launched into a spin tunnel. Such a spin would probably be 
more readily obtainable when the lower rudder is off than when it is on, 
and the ensuing spin would be flatter and faster. Recovery from the 
initial phase of this motion would be satisfactory for either configura- 
tion. However, if the spin is allowed to develop fully, difficulty may 
be experienced in effecting a recovery for the lower-rudder-off config- 
uration, but for the lower-rudder-on configuration satisfactory recovery 
would be obtainable by optimum control technique. 

The results of the investigation also indicated that the magnitude 
of the static pitching moment, the effective dihedral, and the pitch 
damping can be critical in the spin entry and spinning motions for modern 
aircraft similar to the subject configuration and may mean the difference 
between the airplane experiencing a spin or no spin. 
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INTRODUCTION 


For a number of years, the Langley 20-foot free-spinning t unn el 
has been used to test small-scale dynamically ballasted models of air- 
planes. Properly interpreted, the results obtained from tests have 
usually been adequate for predicting the spin and recovery character- 
istics of airplanes represented by the dynamic models. However, modern 
trends in airplane design, primarily very long fuselage forebodies and 
mass distributed very heavily in the fuselage, have caused Reynolds 


number and spin-entry-technique differences between airplanes and models L 

to become extremely important factors which sometimes have marip proper 9 

interpretation of the spin- tunnel-model results quite difficult (ref. l). 5 

Predicting the airplane spin and recovery behavior remains a serious 3 

problem, because modern high-speed airplanes in the categories of 
fighter, interceptor, and research airplanes continue to enter spins 
inadvertently with occasional loss of life and property. 


In an attempt to enable better prediction of full-scale spin and 
recovery characteristics of modern design and to add to the general 
knowledge of how various design and control-sequence factors may affect 
spins and related post-stall transient motions, including incipient 
spins, two additional research and evaluation techniques have been 
initiated at the Langley Research Center to complement the free-spinning- 
tunnel tests. The two new techniques are the use of large radio- 
controlled free-flight dynamic models dropped at altitude and flown into 
various maneuvers conducive of spin entry (ref. 2), and the use of ana- 
lytical techniques which enable calculations of the motions to be made 
on a high-speed automatic digital computer by using static and rotary 
wind-tunnel aerodynamic data and six-degree-of- freedom equations of 
motion (refs. 1 and 3)- 

As regards the analytical technique using a digital computer, inde- 
pendent investigations described in references 3 and 4 for current delta- 
wing and sweptback-wing fighter airplanes, respectively, have indicated 
that it is possible to simulate known airplane spin-entry, developed- spin, 
and recovery motions provided adequate aerodynamic-data coefficients from 
measurements or estimations are used. Developed to its ultimate, the 
analytical technique should enable the prediction of an airplane's spin 
and recovery characteristics before full-scale flights are made. The 
accuracy of such a prediction will, of course, depend on how accurately 
the required aerodynamic data have been obtained in wind tunnels or by 
estimations and how representative are these data of the full-scale 
airplane . 

The analytical study described in this paper was made to aid in the 
prediction of the spin and recovery characteristics of the North American 
X-15 airplane and was conducted in conjunction with free-spinning-tunnel 
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tests made on a l/30-scale dynamic model (unpublished results) and with 
tests on a l/7-scale radio-controlled model dropped from a helicopter 
(ref. 5) . The experimental tests of these X-15 models were made at 
relatively low Reynolds numbers , and it was believed that Reynolds num- 
ber effects on the fuselage forebody might appreciably alter the spin 
and recovery characteristics of this design. 

The approach used in the analytical study was to use low Reynolds 
number aerodynamic force and moment data first to try to simulate the 
spin entry, the ensuing spin, and the recovery therefrom obtained from 
dynamic-model tests. Then, high Reynolds number data from reference 6 
were inserted into the computer and the resulting effects on spin entry 
and spin characteristics were calculated. 


SYMBOLS 


The body system of axes is used. This system of axes, related 
angles, and positive directions of corresponding forces and moments are 
illustrated in figure 1. 

b wing span, ft 

% 

Ct rolling-moment coefficient, 

ipV R %b 

% 

C m pitching-moment coefficient, 

|pVR 2 Sc 

M z 

C n yawing -moment coefficient, 

|fV R 2 Sb 

F x 

C Y longitudinal-force coefficient, 

f y 

C v side-force coefficient, 

|^ r 2 s 
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F Z 

C 7 normal-force coefficient. 

|pV R % 

c mean aerodynamic chord, ft 

F x longitudinal force acting along X body axis, lb 

Fy lateral force acting along Y body axis, lb 

Fg normal force acting along Z body axis, lb 

g acceleration due to gravity, 32.17 ft /sec 2 

altitude at beginning of time increment, ft 
h 2 altitude at end of time increment, ft 

Ix>Iy>Iz moments of inertia about X, Y, and Z body axes, respectively, 

slug-ft 2 

Z maximum depth of fuselage where wing-fuselage intersect , ft 

My rolling moment acting about X body axis, ft-lb 

My pitching moment acting about Y body axis, ft-lb 

Mz yawing moment acting about Z body axis, ft-lb 

m mass of airplane, W/g, slugs 

p,q,r components of resultant angular velocity ft about X, Y, 

and Z body axes, respectively, radians/sec 

R Reynolds number, Vl/v 

S wing area, sq ft 

t time, sec 

u,v,w components of resultant velocity Vp along X, Y, and 

Z body axes, respectively, ft/sec 

vertical component of velocity of airplane center of 
gravity (rate of descent), ft/sec 


V 


hi OMAfO 



resultant linear velocity, ft/sec 
weight, lb 
body axes 

angle of attack, angle between relative wind Vr projected 
into XZ-plane of symmetry and X body axis, positive when 
relative wind comes from below XY body plane, deg 

angle of sideslip, angle between relative wind Vr and 
projection of relative wind on XZ-plane, positive when 
relative wind comes from right of plane of symmetry, deg 

aileron deflection (differential deflection of horizontal 
control) with respect to chord line of wing, positive 
with trailing edge of right aileron down, deg 

horizontal-tail deflection with respect to fuselage refer- 
ence line, positive with trailing edge down, deg 

rudder deflection with respect to fin, positive with 
trailing edge to left, deg 

total angular movement of X body axis from horizontal plane 
measured in vertical plane, positive when airplane nose 
is above horizontal plane, radians or deg 

kinematic viscosity, ft 2 /sec 

air density, slugs /cu ft 

angle between Y body axis and horizontal measured in 

vertical plane, positive for erect spins when right wing 
downward and for inverted spins when left wing downward, 
deg 

total angular movement of Y body axis from horizontal plane 
measured in YZ body plane, positive when clockwise as 
viewed from rear of airplane (if X body axis is vertical, 
0 e is measured from a reference position in horizontal 
plane), radians 

horizontal component of total angular deflection of X body 
axis from reference position in horizontal plane, positive 
when clockwise as viewed from vertically above airplane, 
radians 
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m 


ft resultant angular velocity, radians/sec (or rps where noted) 
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AC^ a incremental rolling -moment coefficient due to aileron 

’ deflection 

AC ^ r incremental rolling -moment coefficient due to rudder 

deflection 

AC nj a incremental yawing-moment coefficient due to aileron 

deflection 

AC n ^ r incremental yawing-moment coefficient due to rudder 

deflection 

ACy^a incremental side-force coefficient due to aileron 

deflection 

ACy^ r incremental side-force coefficient due to rudder deflection 

A dot over a symbol represents a derivative with respect to time; 

for example, u = 

dt 


PROCEDURES AND CALCULATIONS 


Spin entry and spin motions were calculated by a high-speed digital 
computer which solved the equations of motion and associated formulas 
listed in the appendix. The equations of motion are Euler's equations 
representing six degrees of freedom along and about the airplane body 
system of axes. (See fig. 1 for illustration of body axes.) The mass 
and dimensional characteristics used in the calculations are listed in 
table I and a three- view sketch of the airplane Is shown in figure 2. 

Most of the aerodynamic data used were nonlinear and are presented 
in the plots shown in figures 3 to 12. The data in figures 3 to 7 were 
based on measurements obtained in the Langley free-f light tunnel for a 
Reynolds number of about 430,000 based on c - and of about 195,000 based 
on the maximum vertical depth of the model fuselage nose section. The 
data in figures 8 to 10 were obtained In the Langley high-speed 7 _ by 
10-foot tunnel for a Reynolds numb'ef of about 2,000,000 based on c - and 
of about 865,000 based on the maximum vertical depth of the model fuselage 
nose section. Values of the damping-in-pitch derivative C^ were held 

constant ^Cmq = -10.0 as obtained from low-angle-of- attack tests in 

ref. 7 ) for most of the inveptigation. Hovtever, during the investiga- 
tion measured values for the range of a between the stall and 90 ° were 
obtained in a separate wind-tunnel investigation and were used in some 
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of the final calculations of this investigation. These measured values 

of Cjn are shown in figure 12. 

Q. 


as 

is, 

and 


The rotary derivatives presented in figures 7 and 12 were obtained 
combination derivatives which include the effects of £ and a - that 
Cjp is actually sin a^, C nr is actually ^C nr - C n £ cos a 

C T21< ^ is actually /C + C^X. However, inasmuch as the full deriva- 


tives could not be separated into their component parts, it was arbitrarily 
decided for this investigation to treat the derivatives as though they 
were due solely to angular velocities about body axes. These damping 
derivatives, as well as the incremental force and moment coefficients 
due to deflection of the ailerons and rudder, were available only from 
low Reynolds number tests and were arbitrarily assumed to be invariant 
with Reynolds number. 


For the calculations made in an attempt to simulate spin entries 
like that of the radio-controlled model or of an airplane, the effect 
of the variation of air density with altitude was included. For the 
type of calculation simulating conditions under which the model is 
launched into a spin tunnel, air density was maintained constant as it 
is in the Langley 20-foot free- spinning tunnel, and representative 
initial values of attitude and rate of rotation were selected. The 
inputs simulating airplane control movements were introduced into the 
computer by means of appropriate switches. The timing and direction of 
these inputs were based on observed time histories of the computed motion 
as presented by computer print-out tables. 

The significance of motions calculated after application of con- 
trols for attempted recoveries was evaluated on the basis of the fol- 
lowing considerations: An airplane is considered to have recovered 

from a spin when the angle of attack at the center of gravity has gone 
below the angle of attack at which the wing stalls. Usually, when this 
attitude is achieved, the airplane enters a steep dive without rotation 
(r =0). In some cases, however, the airplane may be turning or rolling 
in a spiral glide or an aileron roll. Also, sometimes, the airplane may 
roll or pitch to an inverted attitude from the erect spin and may still 
have some rotation but is out of the original erect spin. 

For convenience, a list of the aerodynamic data used in this inves- 
tigation is presented in table II. 
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BRIEF REVIEW OF DYNAMIC -MODEL-TEST RESULTS 


Results reported in reference 5 and represented for convenience in 
figure 13 indicate that the 1/7-scale dynamic radio-controlled model 
with lower rudder off and its center of gravity at 19.5 percent c was 
flown Into a spin-entry maneuver from trimmed gliding flight and, with 
the differentially operable horizontal tail maintained in an aileron- 
against-the-spin position (stick left in a spin to the pilot's right), 
made a three-turn spin during which the angle of attack varied from 60° 
to 90 ° and in which the rate of spin rotation averaged about 1.1 radians 
per second ( 0.175 nps) for the three turns. Good spin recovery was 
achieved from that spin by reversing the differential tails, hereinafter 
referred to as ailerons, to with the spin, as is also shown in figure 13 . 
Results of free-spinning tunnel tests made at the Langley Research Center 
with a l/30-scale dynamic model with lower rudder off and ailerons against 
the spin indicated that, for center-of -gravity positions ranging from 
0 percent to 20 percent c , fully developed spins were possible in which 
G could be as high as I.85 radians per second (0.335 rps) and in which 
a varied from about 60 ° to 90 ° • From these spins, recoveries varied 
from fast to slow and, on this basis, recovery characteristics were con- 
sidered unsatisfactory (lower rudder off). 


RESULTS AND DISCUSSION 


Calculated results of the analytical study are presented in fig- 
ures 14 to 26 as time histories of angle of attack a, angle of wing 
tilt 0, yawing velocity r or resultant angular velocity fl, control- 
surface positions, and spinning turns completed. Although only these 
few pertinent variables of the motions are presented, time histories 
of all the attitude, velocity, and acceleration variables In the equa- 
tions of motion were obtained. 


Low Reynolds Number Calculations With Lower Rudder Off 

Calculations were made to simulate the result (fig. 13 ) which was 
obtained during a flight of the 1/7- scale radio-controlled model at low 
Reynolds number and with the lower rudder off. Zero-time conditions 
for the calculations were estimated from the results of the model test 
and are listed in column A of table III. The initial calculated results, 
shown in figure l 4 , did not agree with the model results (fig. 13), but 
rather indicated a continuing cyclic motion in pitch with variations of 
angle of attack between 20° and 80° and with a rate of rotation averaging 
about 0.7 radian per second. Because of the large oscillations in pitch, 
it was reasoned that perhaps some adjustment was desirable to the 
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magnitudes of aerodynamic pitching moments used. From an examination of 
some pitching-moment data plotted against sideslip angle for this con- 
figuration (ref. 5)> it was apparent that large reductions in negative 
values of C m could occur when large sideslip angles are present, such 
as those obtained at high angles of attack during the calculated attempted 
spin entry. Considering this fact and also the assumption that shielding 
of the horizontal tail by the vertical stub and fuselage area under the 
horizontal tail would possibly reduce the aerodynamic nose-down pitching 
moment in a spinning motion, it was decided to investigate the effect of 
reducing the curve of against a as an input into the equations of 

motion. Arbitrarily, the static-pitching -moment values (fig. 3) were 
reduced by approximately 50 percent and another calculation was made. 

The results are shown in figure 15, and after three turns of this cal- 
culated spin motion the angle of attack and rate of rotation agreed with 
those for the three- turn spin obtained during the model tests (fig. 13). 

A recovery from this spin (fig. 15) was attempted by reversing the ailer- 
ons and rudder, and recovery was obtained in approximately one turn, a 
result which also is in qualitative agreement with the experimental 
result . 

It may be noted that included among the initial conditions for the 
spin-entry calculation just discussed (initial conditions are given in 
col. A of table III) was a nose-down attitude, 0 e = -60°, used to simu- 
late the experimental model initial attitude. As a matter of interest, 
a calculation was made in which an attempt to enter a spin was started 
from an arbitrary nose-up attitude, 0 e = 20°. (See table III, col. B.) 
The results are shown In figure l6 and indicate a spin that for the first 
three turns was also generally similar to the experimentally obtained 
spin, as had been the three -turn spin calculated starting from the 
initial nose-low attitude. The results from both these calculated three- 
turn spins are considered to have simulated the experimental 1/7-scale- 
model spin adequately. 

As may be seen in figure 16, no recovery was attempted from the 
second of the two calculated three- turn spins. Instead, pro- spin con- 
trols were maintained and, as can be seen, the spin-rotation rate con- 
tinued to increase until after a total of eight spinning turns. At that 
time, Q was averaging about 2.1 radians per second (0.335 rps). Per- 
haps if it had been practicable to allow the experimental model to spin 
longer, it too may have eventually increased its rate of rotation. At 
any rate, inasmuch as the l/30-scale free- spinning-tunnel model had 
exhibited developed spins rotating as high as 1.85 radians per second 
(0.295 rps), as noted earlier, the calculated and experimental (spin 
tunnel) developed spins are considered to be in qualitative agreement, 
even though the calculated value of Q was somewhat higher than the 
experimentally obtained value of Q. As regards recovery character- 
istics from the developed spin, a recovery attempt by reversing ailerons 
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and rudder was made from the calculated spin after eight turns (fig. l6) 
and the result indicated no recovery in four turns. This is considered 
to he in qualitative agreement with the free- spinning- tunnel experimental 
result, which (as noted before) indicated unsatisfactory recovery char- 
acteristics with lower rudder off because recoveries varied from fast to 
slow. 


High Reynolds Number Calculations 

Inasmuch as it appeared at this point in the analytical study that 
the experimental model motions had been reasonably simulated, the low 
Reynolds number aerodynamic data were removed from the computer and cor- 
responding high Reynolds number inputs were inserted for the attempts 
to calculate possible airplane spin characteristics. The aerodynamic 
data that were changed because of Reynolds number differences appear in 
figures 8 and 9 as plots of C m , C z , C x , C Yft j C n , and Cj against 

a. As mentioned previously, the damping derivatives and the incremental 
force and moment coefficients due to deflection of ailerons and rudder, 
all measured only at low Reynolds number, were not changed for the high 
Reynolds number calculations. Also, because the aerodynamic -data meas- 
urements indicated no significant difference in the damping coefficients 
between the lower-rudder-on and lower-rudder-off conf igurations in the 
high angle-of -attack range, the values used in the calculations for 
these damping coefficients are the same for lower rudder on and for 
lower rudder off. Incremental effects of the lower rudder on the static 
derivatives C n ^ and C^, as v ell as on the incremental moments due 

to rudder deflection, were appreciable, however, and these effects were 
included in the calculation inputs. 

Lower rudder on, simulated launching with rotary motion .- In order 
to determine whether the configuration with lower rudder on would spin 
at high Reynolds number, a launching with spin rotation was simulated. 

The Initial conditions used for this launching were obtained from free- 
spinning-tunnel model tests (table III, col. C). The calculated motion 
obtained is presented in figure 17 and indicates Increasing oscillations 
in both pitch and roll as the rate of spin rotation increased. The 
motion oscillated out of range of the computer’s ability to continue 
the calculation and no developed spin was achieved. Inasmuch as it had 
been found in the low Reynolds number calculations that 100 percent of 
the static pitching moment probably does not act during spins, the curve 
of C m against a (fig. 8) was reduced approximately 50 percent and 
another calculation was made. The resulting motion is in figure l8 and, 
as can be seen, a very steady spin is indicated. In order to see whether 
a lesser reduction in the measured static pitching moment would allow 
the attainment of a spin, other calculations were made when 85 percent 
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and 70 percent of the measured static values of „vere used. These 
calculations shoved that with An = 85 percent there was no spin 
obtained, vhereas with C m = 70 percent there was a spin. As shown in 
figure 19, where An values approximately equal to 70 percent of the 


measured static values were used, the rate of rotation is higher than 
that obtained with An values of 50 percent (fig, l8), and such a 
result can be expected if a spin can be maintained (ref* l). Also, a 
recovery was attempted from this spin motion (fig. 19) , and it can be 


seen that the airplane recovered in approximately 2i turns after antispin 


controls were deployed (rudder and ailerons reversed). 


In view of the aforementioned results, it might be said that the 
full-scale airplane with the lower rudder on would spin if It was 
launched into a spin and if the pitching moment acting during spins was 
no more than about 70 percent of the measured static values shown in 
figure 8. 


Lower rudder on, simulated spin-entry motion .- Inasmuch as a spin 
was obtained after a launching with rotation, as indicated for the high 
Reynolds number case when values equal to 70 percent of the measured 
static values of An were used, some attempts were made to obtain a 
spin from a more normal entry maneuver for the configuration with lower 
rudder on. For this calculation, the initial conditions shown In 
table III, column D, were used together with the same set of aerodynamic- 
data inputs as was utilized for the spin- launching calculation when 
An ~ 70 percent of the measured static values (fig. 19)* The resulting 
motion is shown in figure 20 and indicated a diverging oscillation which 
ended in the airplane experiencing no spin as the angle of attack went 
below zero. Although not presented, a calculation was also made with 
An ~ 50 percent of the measured static values and the results obtained 
were similar to those obtained with C m values of 70 percent. 

In an attempt to understand better the results just discussed, 
another spin-entry calculation was made in which larger negative values 
of Cjp were used. (See fig. 4 for values.) Seventy percent An 

values were also used. The results are shown In figure 21 and indicate 
that a spin was now obtained from entry. The angle of attack during this 
spin was approximately 50° and r was about 0.7 radian/second. From 
these results, it appears that the X-15 airplane with lower rudder on 
will not enter a spin from this type of entry because of the magnitude 
of the dihedral effect. 


Lower rudder off, simulated launching with rotary motion .- In order 
to investigate the spin characteristics, at high Reynolds number, of the 
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lover -rudder -off configuration, a calculation vas made in which the air- 
plane was launched with applied spin rotation and with (^^70 percent 

of the measured static values. The result is shown in figure 22. A 
comparison of figures 22 and 19 indicates that the lower-rudder- off spin 
is flatter and faster than the lower-rudder-on spin, as would be expected. 
It can also be seen that the airplane did not recover from this fast- 

rotating spin (fig. 22), whereas it had recovered In 2^ turns for the 

2 

slower rotating spin shown in figure 19 for the lower-rudder -on case. 

Lower rudder off, simulated spin-entry motion . - In an attempt to 
obtain a spin from a more normal type of entry for the lover-rudder-off 
configuration by using high Reynolds number data, the initial conditions 
shown in table III, column D, and values ’JO percent of the measured 

static values were used. The results indicated that the airplane experi- 
enced no spin in that the angle of attack went below zero. Another sim- 
ilar calculation was made with « 50 percent of the measured static 
values and again no spin was indicated. 

At this point a calculation was made to see what effect further 
reduction in the aerodynamic -pitching-moment input would have on the 
tendency to enter a spin. The values shown in figure 11 were used 

and the results, shown In figure 25, indicate that a spin was attained. 
Although not investigated, it is believed that, if larger negative values 
with Cjn values equal to 70 percent of the measured static 

values had been used, as had been done for the lower-rudder- on calcula- 
tions, they would have also enabled the attainment of a spin. In view 
of these considerations and some additional unpresented calculation 
results, indications are thal} the magnitudes of static pitching moment 
and effective dihedral have important interrelated effects as regards 
the determination of whether the subject configuration will spin. It 
appears that, when relatively large nose-down aerodynamic pitching 
moments are present, corresponding large amounts of effective dihedral 
are required for the attainment and maintenance of spins. 

From the results just discussed, it appears that even the lower- 
rudder-off configuration at high Reynolds number will not enter a spin 
from gliding- flight attitudes. Factors contributing to the prevention 
of spins would be the magnitudes of the pitching and rolling moments 
present. However, if the airplane could, for some reason, achieve an 
attitude and motion simulating a rotary launching, such as Is obtained 
in a free -spinning tunnel, the airplane could probably enter a spin. 
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Importance of Magnitude of Pitch Damping 

As stated previously, a constant value of -10.0 (based on results 
in ref. 7) was used for C m ^ in most of the calculations. It was 

decided to investigate briefly the effects of varying C^ on the spin. 

The spin motion shown in figure 22 was arbitrarily chosen. With the 
same conditions used as inputs, but with = ~5»0, a calculation was 

made and the resulting motion is shown in figure 24. A comparison of 
figures 22 and 24 indicates that this reduction in C^ had little or 

no effect on the spinning motion when the airplane was launched with 
applied spin rotation. 

In order to see whether a reduction in Cmq would have any effect 

on a spin motion obtained from entry without applied rotation, calcula- 
tions were made which had the same initial inputs as shown in figure 23 
except for different values of Cni . The first arbitrary value used was 

C~ = -7.0. The results are shown in figure 25 and indicate that now a 

m q 

spin could not be attained. Because of the indicated importance of the 
magnitude of C mq as affecting spin entry, some wind-tunnel tests were 

made (as previously indicated) to determine the variation of C^ with 

angle of attack for this design. The measured results are shown in fig- 
ure 12. A spin-entry calculation, without applied rotation, was made by 
using this curve as an input; the results are given in figure 2 6 and 
indicated that a spin motion was again attained and that this spin was 
similar to the one In which the constant Cm^ = -10.0 was used (fig. 23). 

It may be concluded from these results that the values of used as 

an Input can make a difference between the airplane experiencing a spin 
or no spin. Indications are that carefully measured values of this 
derivative for angles of attack ranging from 0° to 90 ° are desirable in 
analytical spin studies. 


CONCLUDING REMARKS 


The North American X-15 airplane with the lower rudder on or off is 
not likely to enter a spin from trimmed gliding flight. A developed spin 
might be obtained if the airplane receives some violent disturbance which 
would tend to put it in a flight condition having rotation and a high 
angle of attack, approximating the attitude and motion of the model when 
launched into a spin tunnel. Such a spin would probably be more readily 
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obtainable when the lower rudder is off them when it is on, and the 
ensuing spin would be flatter and faster. Recovery from the initial 
phase of this motion would be satisfactory for either configuration. 
However, if the spin is allowed to develop fully, difficulty may be 
experienced in effecting a recovery for the lower-rudder-off configura- 
tion, but for the lower-rudder-on configuration satisfactory recovery 
would be obtainable by optimum control technique. 

The magnitude of static pitching moment, effective dihedral, and 
pitch damping cam be critical in the spin entry and spinning motions 
for modern aircraft similar to the subject configuration and may mean 
the difference between the airplane experiencing a spin or no spin. 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va., March 1 6, i960. 
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APPENDIX 


EQUATIONS OF MOTION AND ASSOCIATED FORMULAS 


The equations of motion used in calculating the spinning motions 


were: 


. Iy * 1 2 pV R Sb pV R Sb pV R Sb pV R Sb 

p = — qr + — — C 7 „3 + — C 7 p + — ^ AC 7 ^ + AC 


— qr + + °y + 


2I\ 


l, a. 


I Z " I X 


2o— 


pV R -Sc P V R Sc* 


q pr + c m + jjV" c mq 


2I\ 


I x - I y pV R Sb pV R Sb 2 pV R 2 Sb 


r = 


Iz 


pq + 


2IZ 


C np P 


kl'? 


Cnr 1 * 


ACn . t* + 


pV R ^Sb 


ot "n,r 

2I Z 2I Z 


AC 


n,a 


u = -g sin 0 e + vr - wq + 


pV R 2s 


2m 


P V R 2 S P V r 2 S P V r 2 S 

V = g cos 0 e sin 0 „ + wp - ur H Cy D 0 + ACy >. + ACy , 

e re * 2m *3 2m x > r 2m lfl 


, pV R 2s 

w = g cos 0 e cos p e + uq - vp + C7 

2m 


In addition, the following formulas were used: 


a = tan"^- — 
u 


V-M vji vo tr 1 
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P = sin" 1 ^ 

V R 

V R = \j u 2 + v 2 + v 2 

« = ^p 2 + q 2 + r 2 

V = -u sin 0 e + v cos 0 e sin 0 e + w cos 0 e cos 0 e 

h2 = h]_ - AtV 
0 e = q cos 0 e - r sin 0e 
0 e = p + r tan 0 e cos 0 e + q tan 0 e sin 0 e 

#e “ p 
xjr = 

e sin 0 e 

f*e dt 

Turns in spin = 

dA 

j . -1 sin 0 

0 = sin x 11 

cos 0g 
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TABLE I.- MASS AND DIMENSIONAL CHARACTERISTICS 


Values given in terms of full-scale airplane 


High Reynolds number 
calculations (from 
airplane data) 

10.27 
22.36 
200 
13,577 

20.0 

3,823 

75,313 
76,738 

Maximum control deflections: 


S r , deg ±7-5 ±7-5 

5 a > de 6 ±7-5 ±7-5 


Low Reynolds number 
calculations (same 
as l/7-scale 
dynamic model) 


C , ft 10.27 

b, ft 22.36 

S, sq ft 200 

W, lb 12,575 

Center-of -gravity loca- 
tion, percent c . . . 19*5 

I}£, slug-ft 2 4,288 

l Y , slug-ft 2 73,384 

I z , slug-ft 2 ..... 74,867 
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TABLE II— LIST OF AERODYNAMIC DATA 


Coefficient 

Reynolds 

number 

Basic configuration with - 

Figure 

" C Z 

against a 

Low 

Lower rudder off j 

c.g. = 19.5 percent cj 
6 h = -30° 

3 

Cl V c v c Yp 

Low 

Lower rudder off 

k 

against a 




AC n ^ a> ^Y,a 
against a 

Low 

Lower rudder on and offj 
& a = ±7|° 

5 

^n,r> ^Y,r 
against a 

Low 

Lower rudder on and offj 

S r = ±t| 

6 

C ip and C nr 

Low 

Lower rudder off 

7 

against a 




c m> " C X> ~ C Z 
against a 

High 

Lower rudder onj 

c.g. = 20 percent 5; 

5 h = -35° 

8 

Cn p> c v % 

High 

Lower rudder on 

9 

against a 




% <“ d % 
against a 

Cjn 1 against a 

High 

Lower rudder off 

10 

11 



C^ against a 

335,000 

Lower rudder offj 6^ = -35° 

12 


■^Arbitrary values of Cm used. 


h-3 0\ IT\ K"\ 










vx Ui vo tr 1 


21 


TABLE III.- INITIAL CONDITIONS USED IN CALCULATIONS 



Column A 

Column B 

Column C 

Column D 

<x, deg 

15 

20 

80 

15 

P, cleg 

-15 

-15 

-1 

-15 

0 e , deg 

-6o 

.20 

-10 

-6o 

0/ deg 

-5 

-9 

0 

-5 

de e 

0 

0 

7-5 

0 

5 a , deg 

0 

0 

- 7-5 

0 

u, ft/sec 

246 

240 

445 

246 

v, ft/sec 

-68 

-68 

-4 

-68 

v, ft/sec 

66 

8? 

251 

66 

p, radians/sec 

0 

0 

.245 

0 

q, radians/sec 

.7 

• 3 

0 

• 7 

r, radians /sec 

0 

0 

1.379 

0 

hi, ft 

19,000 

19,000 

50,000 

50,000 



Horizontal 



(c) 9 e and = 0, and in this case 0 = 0 e 


Figure 1.- Body system of axes and related angles 









40 50 

of attack, a, deg 


60 


ning-moment, axial-force } 
attack. Center of gravj 







Angle of attack, a, deg 

Figure 6.- Variation in increments in the lateral force and moment 
coefficients with angle of attack due to deflecting rudder. 















Time , sec (ful I scale) 


Figure 13.- Results of radio-controlled dynamic -model test for lower 
rudder-off configuration. (Data from ref. 5-) 
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Figure 14 . - Results of initial calculated spin-entry attempt. Low 
Reynolds number; lower rudder off. 
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Figure 15.- Spin entry, spin, and recovery calculated for low Reynolds 
number. Lower rudder off; C m approximately 50 percent of measured 

static values. 
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Figure 17.- Calculation simulating spin- tunnel model prerotated launching 
technique at high Reynolds number. Lower rudder on. 
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Figure 20.- Calculated spin-entry attempt at high Reynolds number. 

Lower rudder on; C m approximately 70 percent of measured stati 
values . 
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Figure 21 . - ~ Calculated spin entry and spin for high Reynolds number 
Lower rudder on; C m approximately JO percent of measured static 
values; larger negative C 7 values used (see fig. 4). 
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Time, sec (full scale) 


Figure 25*- Calculated spin-entry attempt at high Reynolds number. 
Lower rudder off; arbitrarily reduced C m values used (see 

fig. 11); = -7.0. 
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Figure 26.- Calculated spin entry and spin for high Reynolds number. 

Lower rudder off; arbitrarily reduced C m values used (see fig. 11); 
curve of C m ^ against a in figure 12 used. 
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